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LONG-TERM  GOALS 

The  long-term  goal  of  the  REVEAL  project  is  to  develop  and  demonstrate  a  passive  sonar 
signal  processing  structure  that  exploits  available  knowledge  of  the  environment,  including 
uncertainty,  to  detect,  localize  and  classify  targets. 

OBJECTIVES 

The  FY12  objectives  were  to: 

•  Continue  development  and  evaluation  of  a  striation-based  processor  for  passive  sonar 
localization,  and  apply  the  algorithm  to  the  south  Florida  array  CALOPS  data. 

•  Apply  the  mode  space  processor  developed  by  Premus  to  estimate  passive  source  depth  using 
horizontal  and  vertical  line  arrays  to  data  from  the  NURC  REP1 1  B  data,  and  evaluate  sensitivity 
to  key  environmental  factors. 

APPROACH 

The  primary  effort  of  this  project  has  been  to  develop  and  evaluate  passive  sonar  receiver 
structures  that  make  use  of  knowledge  of  the  environment  and  received  signal  characteristics  to 
detect,  localize  and  classify  the  passive  source,  particularly  in  shallow  water.  Several  structures 
have  been  or  are  being  investigated. 

In  shallow  water,  passive  sonar  context,  the  characteristics  of  received  signals  are 
determined  by  the  properties  of  the  source  and  by  propagation  through  the  ocean.  Our 
hypothesis  is  that  knowledge  of  the  environment  can  be  used  to  estimate  propagation  effects, 
which  then  provide  information  about  the  source  location.  Because  the  ocean  environment  is 
dynamic  and  variable  in  time  and  space,  a  statistical  approach  is  necessary. 

WORK  COMPLETED 

In  a  shallow  water  waveguide,  where  the  distance  between  an  acoustical  point  source  and 
receiver  is  several  times  greater  than  the  water  depth,  the  resulting  pressure  field  may  be 
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expressed  by  a  finite  sum  of  normal  modes.  A  superposition  of  modes  represents  the  sound 
pressure  at  the  receiver,  i.e. 

N 

P(r,zr,zs)  =  YJ  Bn  (zr ,  zs )  exp  [ikrn  r]  , 

n= 1 


(l) 


when  the  depth  of  the  source,  Z5 ,  and  range,  r  ,  from  the  omnidirectional  point  receiver  located  at  depth. 


zr,  are  known  and 


Bn(-r^s)  = 


k  2  TV  ^ 


1/2 


krn>-J 


(2) 


represents  the  contribution  from  each  mode.  The  mode  eigenfunctions  and  eigenvalues  are, 

|  //  J  ^ ,  k 

YnK  '  and  rn ,  respectively,  and  rn  is  the  horizontal  component  of  the  wavenumber  for  mode  m. 
Lastly,  the  number  of  modes  contributing  to  the  pressure  field  ™  is  a  function  of  source 
frequency  and  waveguide  depth. 

If  the  acoustical  source  frequency  is  high  enough  to  excite  more  than  one  propagating  mode 
in  the  waveguide,  the  modes  will  interact  with  one  another  to  create  an  interference  pattern  in 
range,  depth,  and  frequency.  Writing  the  difference  in  mode  wavenumbers  as 
A  knm((6)  =  krm  (a>)  -  knl(co ) ,  where  m*n,  the  received  pressure  spectrum  can  be  written  as 


I(r,zr,zs,co)  =  \P\2  =YjYjBn,Bn  exp[/AA'nw(ru)r].  (3) 

m  n 

Equation  (3)  describes  the  intensity  interference  pattern  between  different  modes  excited  by  a 
broadband  noise  source.  This  interference  pattern  can  be  seen  as  spectral  striations,  which  have 
unique  character  determined  by  the  relative  location  of  the  source  and  receiver  in  the  waveguide, 
as  well  as  the  parameters  that  describe  the  waveguide. 

Equation  (3)  also  provides  a  method  of  determining  the  full  intensity  field  based  upon  the 
waveguide  parameters,  however,  it  is  useful  to  relate  these  parameters  to  lines  of  constant 
intensity  within  the  full  field.  Setting  the  derivative  of  Eqn.  (3)  equal  to  zero,  a  constant  intensity 
contour  may  be  described  as 

r  a  r 

dI(r,co)  = - A  ftH - Ar  =  0  (4) 

dco  dr 

where  A  0)  —  o>  —  0)(]  and  A r  =  r  —  r0.  Using  Eqns.  (3)  and  (4),  and  solving  for  the  slope  of  the  contour  in 
the  range-frequency  plane  yields 


A eo  _  ^knJo))l(o  =  o)(]  &S""(o)) 

A r  r0  dAknin((o)/d(o  r0  AS™"(co)  ' 

The  contours  are  typically  line-like  and  are  often  referred  to  as  striations.  In  Eqn.  (5),  the  mode- 
dependent  wavenumber  terms  have  been  identified  as  group  and  phase  slowness,  Sg  and  So .  Using  the 
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WKB  approximation4,  it  is  possible  to  show  that  the  relationship  between  S„  and  Sp  does  not  depend  on 
either  mode  number  or  frequency.  This  leads  to  identification  of  the  waveguide  invariant,  /? ,  defined  by 
Chuprov1  as 


(6) 


The  waveguide  invariant  is  a  useful  parameter  for  range  estimation  because  it  relates 
striation  slope  taken  from  observed  interference  patterns  to  r0 ,  the  range  between  the  source  and 

receiver  and  co0 ,  the  center  frequency  of  the  band  in  which  the  striation  is  observed.  More 
detailed  derivations  of  the  waveguide  invariant  may  be  found  in  other  references I’2,5,6,7. 

RESULTS 


The  algorithm  just  described  was  applied  to  data  from  an  experiment  which  took  place  along 
the  continental  shelf,  east  of  Port  of  the  Everglades,  Florida.  Acoustical  data  were  taken  using  a 
1 1 8  element  horizontal  line  array  (HLA)  positioned  on  the  sea  floor  in  approximately  260  m  of 
water.  The  HLA  was  approximately  220  m  long  and  was  oriented  from  west  to  east.  In  addition, 
the  R/V Seward  Johnson  transited  northward  from  the  array  with  an  acoustical  source  towed  at  a 
depth  of  1 00  m.  Although  the  towed  source  is  not  of  interest  in  this  paper,  the  Seward  Johnson, 
recorded  data  using  her  radar's  Automatic  Identification  System  (AIS),  which  provided  position 
data  for  surface  vessels  within  25  km  of  the  HLA.  The  AIS  ship  position  data  were  used  as 
ground  truth  in  the  localization  tests  presented  below. 

Over  the  course  of  the  experiment,  numerous  sound  velocity  profile  (SVP)  measurements 
were  made  from  the  Seward  Johnson  using  stationary  casts  of  a  conductivity,  temperature,  and 
depth  (CTD)  sensor  and  moving  casts  with  expendable  bathythermographs  (XBT).  The 
sediment  parameters  used  for  modeling  wrere  taken  from  Greene  et  al.3  Although  sediment 
parameter  data  exist  for  this  area,  they  are  by  no  means  complete  and  thus  add  uncertainty  to  the 
results  presented  in  the  next  section. 

In  order  to  estimate  range  using  spectral  striations  seen  in  range- frequency  surfaces,  Eqn.  (6) 
must  be  solved  for  source  range,  r0 ,  to  yield 


where  the  striation  slope  is  £  =  ■  The  source  is  at  angle  (f)  from  HLA  broadside  direction  and  is 

sufficiently  far  from  the  array  such  that  received  signals  may  be  approximated  as  plane  waves.  The  angle 
to  the  source  may  be  determined  using  conventional  beamforming.  Spectral  processing  then  is  done  for 
each  channel  in  the  HLA  in  order  to  obtain  the  range-frequency  interference  pattern  along  the  array. 
Striation  slope  may  be  discerned  from  the  spectrograms,  as  shown  for  example  in  Fig.  (1).  This  method, 
however,  is  limited  in  cases  where  the  source  azimuth  angle  is  close  to  broadside,  otherwise  the  projected 
aperture  may  be  insufficient  to  observe  variations  in  frequency  spectra  along  the  array. 

The  remaining  part  of  Eqn.  (6)  that  must  be  estimated  is  the  waveguide  invariant,  /? .  For 
cases  where  the  shallow  water  environment  is  poorly  known,  and  the  source  is  known  to  be  near 


4 


the  surface,  a  value  of  one  can  be  used  for  / 31.  However,  this  is  a  poor  approximation  for  some 
scenarios  because  the  SVP,  receiver  depth,  and  bathymetry  between  the  source  and  receiver 
affect  the  value  of  the  waveguide  invariant.  Any  one  of  a  number  of  Normal  Mode 
methods1"  '5  6  may  be  used  to  determine  /3  depending  on  the  type  of  environment  and  the  degree 
to  which  its  features  are  known. 
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Example  Bilateral  Filtered  Spectrogram 


Figure  1:  Spectrogram  showing  striations  in  the  receive  intensity. 

During  the  eighth  hour  of  acoustical  "Run  1  North"  in  the  CALOPS  experiment,  the  cargo 
ship  Socol  2  passed  just  west  of  the  HLA  traveling  northward.  Her  closest  point  of  approach 
(CPA)  to  the  array  was  approximately  2.5  km,  and  a  sufficient  amount  of  broadband,  noise  was 
radiated  to  excite  the  waveguide  and  create  visible  spectral  striations  along  the  array.  Figure  (1 ) 
shows  the  character  of  the  interference  patterns  seen  during  this  part  of  the  dataset.  These 
striations  are  quite  clear,  however,  current  effort  is  focussed  on  producting  interference  patterns 
which  are  clearly  visible  and  useable  for  slope  estimation. 

In  order  to  accurately  resolve  the  striations,  a  number  of  image  processing  schemes  were 
investigated.  The  spectrograms  were  first  normalized  at  each  frequency  to  remove  completely 
horizontal  features.  They  were  then  processed  using  a  bilateral  filter*  to  accentuate  larger 
features  (e.g.  striations)  and  remove  speckle.  The  result  was  spectrograms  with  far  more 
apparent  striations,  when  striations  were  present.  Even  with  the  aforementioned  processing,  the 
resulting  spectrograms  showed  substantial  fading  over  the  length  of  the  HLA.  The  final 
processing  step  was  to  average  10,  1  second  snapshots  to  create  the  an  image  used  for  slope 
estimation. 
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Estimated  striation  slopes  were  used  in  Eqn.  (6)  to  yield  the  range  estimates  shown  in  in 
Fig.  (2).  In  addition  to  showing  range  estimates  using  P  =  1 ,  the  Range  Dependent  Waveguide 
Invariant  Distribution  (RaDWID)7  method  was  used  to  estimate  p  based  on  known 
environmental  parameters.  Using  RaDWID,  estimated  values  for  p  were  closer  to  0.8  and  thus 
decreased  the  estimated  range  values.  At  close  range  the  range  estimates  are  within  +/-  500  m  of 
true  range.  However,  as  the  source  moves  away  from  the  HLA  and  the  actual  striation  slope 
decreases,  small  errors  in  slope  estimation  result  in  larger  range  estimation  errors.  Also  of  note 
are  the  gaps  present  before  minute  36,  between  minutes  42  and  46,  and  after  minute  52.  These 
were  instances  where  striations  were  not  apparent  in  the  spectrograms. 


SOCOL  2  Acoustical  Range  Estimate  and  AIS  Ground  Truth 


Figure  2:  Estimated  source  range  (circles)  and  ground  truth  (line). 

We  discuss  two  possible  mechanisms  which  might  explain  the  absense  of  visible  striations 
during  certain  intervals.  First,  before  minute  36,  the  ship  aspect  angle  was  just  off  of  the  bow 
and  well  forward  of  the  beam.  Thus  the  hull  of  the  ship  might  block  propulsion  noise  (i.e.  from 
the  propeller  and  machinery).  Second,  the  range  from  the  ship  to  the  HLA  was  lower  between 
minute  42  and  48, and  it  is  possible  that  the  direct  path  was  much  stronger  than  the  multipath 
signal,  resulting  in  only  weak  interference. 

IMPACT/APPLICATIONS 


6 


This  research  appears  to  offer  a  new  automation  capability  for  localizing  passive  sonar  sources.  The 

algorithms  can  be  applied  to  any  undersea  sensor,  such  as  a  submarine  hull  or  towed  array  or  a  fixed 

sensor  mounted  on  the  ocean  floor. 

TRANSITIONS 

The  striation-based  localization  algorithm  was  been  proposed  for  the  FY  13  submarine  Advanced 

Processing  Build  (APB  13)  program  but  was  not  funded.  The  algorithm  will  be  proposed  again  for 

APB15. 

RELATED  PROJECTS 

None. 
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